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Colloidal particles grafted with single-stranded DNA (ssDNA) chains can self-assemble into a number of differ¬ 
ent crystalline structures, where hybridization of the ssDNA chains creates links between colloids stabilizing 
their structure. Depending on the geometry and the size of the particles, the grafting density of the ssDNA 
chains, and the length and choice of DNA sequences, a number of different crystalline structures can be fab¬ 
ricated. However, understanding how these factors contribute synergistically to the self-assembly process of 
DNA-functionalized nano- or micro-sized particles remains an intensive field of research. Moreover, the fabrica¬ 
tion of long-range structures due to kinetic bottlenecks in the self-assembly are additional challenges. Here, we 
discuss the most recent advances from theory and experiment with particular focus put on recent simulation 
studies. 

Keywords: DNA-functionalized nano-particles, self-assembly, experiment, theory, computer simulation 
PACS: 82.70.Dd, 87.U.gk, 81.16.Dn, 61.46.Df, 61.50.Ah, 75.75.Fk 


1. Introduction and overview 

The miniaturization of technology demands new highly efficient methods to fabricate materials with 
desired optical, electronic and mechanical properties. To this end, the self-assembly of nano- and micron¬ 
sized particles functionalized with DNA strands has attracted much attention as a promising candidate 
for the formation of nano-crystals with long-range periodicity (iHSl. as well as in applications such as 
sensing, imaging and therapeutics The role of DNA strands attached onto the surface of the par¬ 

ticles is to stabilize soft non-compact crystalline structures by acting as links between the particles or 
as the intermediate matter. The structures depend on the size and geometry of the nano-particles (NPs), 
the surface density of DNA strands grafted onto the particles and parameters related to the DNA strands, 
such as the number of complementary bases which are available for hybridization (T3|. Additionally, the 
occurrence of kinetically trapped structures during the self-assembly depends on thermodynamic condi¬ 
tions. 

The idea of using DNA strands to create non-compact crystalline nano-structures was introduced by 
Alivisatos et al. G1 and Mirkin et al. GU, where gold (Au) NPs were used. These NPs are also referred 
to as polyvalent DNA Au NP conjugates (16) due to the large number of oligonucleotide chains grafted 
onto the NPs. Since then, much research is dedicated to understanding the self-assembly process of DNA- 
functionalized NPs and designing strategies that lead to specific crystalline structures. After almost two 
decades of intensive research in this field, a very good understanding has been achieved through a num¬ 
ber of most interesting results. However, significant problems remain unsolved; for example, kinetic 
bottlenecks during the self-assembly process may result in the formation of unwanted crystalline, or 
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disordered structures. Moreover, synthesis of DNA-coated NPs in a controlled and accurate manner re¬ 
quires further improvement (13. Therefore, the manufacturing of long-range ordered nano-structures 
required for many applications remains challenging. Additionally, bulk non-compact materials with a 
photonic band-gap in the range of visible wavelengths have not been achieved either (Ts). 

Experimental, theoretical and computer simulation methods have contributed substantially to this 
area of research. While computer simulations offer some advantages related to the control of conditions 
over theory and experiment, limitations in system size and availability of reliable force fields underline 
the importance of a combined approach in studying the self-assembly of DNA-functionalized NPs. This 
review attempts to briefly present the key results in this research area with an emphasis put on the 
simulation side. 

Experimental research and a range of applications of DNA-functionalized NPs have been previously 
discussed in a number of reviews (THTOl ri9lf25l . while very little discussion has been reported on the 
recent progress of theoretical and simulation contribution to this area of research I26H29I . Part of the ex¬ 
perimental discussion has focused on the synthesis, characterization and phase behavior of DNA-coated 
NPs assemblies (lO] [iTl IT9H221 [24l [23, while others have concentrated on the possible applications of 
DNA-functionalized colloids (THU- 

Storhoff and Mirkin (23 have described the synthesis of DNA-programmed materials (30) referring 
to meso- and macroscopic organic structures from DNA at its infancy. Later, Macfarlane et al. have syn- 
opsized the key findings of their work suggesting fundamental design rules to predict the synthesis of 
distinct colloidal crystal structures with lattice parameters within the range 2S-1S0 nm based on the 
independent adjustment of particle size, lattice periodicity and inter-particle distance (T3|23|3ll- For 
example, it was suggested (23 |3T) that the overall hydrodynamic radius of a DNA NP dictates the as¬ 
sembly and packing behavior of the system, where the hydrodynamic radius is the size of the NP with 
the oligonucleotide components. NPs of equal hydrodynamic radius with self-complementary sequences 
self-assemble into fee lattice, while CsCl or bcc lattices are obtained by using two different NPs with com¬ 
plementary DNA sequences (T3, where the hydrodynamic radius can be measured by techniques, such as 
fast Brownian motion analysis (33- It has been found that DNA NPs with equal hydrodynamic radii will 
maximize the number of hybridization events between neighboring colloids. Additionally, the size ratio 
(expressed through the ratio of the hydrodynamic radii) and the DNA hnker ratio (the ratio of the number 
of DNA linkers) between two NPs control the thermodynamically favored crystal structure (23- In partic¬ 
ular, the size ratio determines the packing of particles in the crystal and its stability, while the linker ratio 
determines the stability of the crystal. Therefore, two systems with the same size and linker ratio exhibit 
the same “thermodynamic product” (T3[23- The final crystal structure is the one that maximizes the 
DNA sequence-specific hybridization interactions. It has been also shown that similar structures can be 
produced by varying the rate at which individual DNA linkers dehybridize and subsequently hybridize 
(23l . These conclusions have been reproduced by using a simple model based on the assumption that 
DNA sticky ends should physically contact each other to hybridize and that any sticky ends coming into 
contact will eventually form a DNA duplex. Hence, this simple model can be used as a tool to inform 
the self-assembly of DNA-functionalized spherical NPs (23- The work of Macfarlane et al. (TTl 1231131) 
highlights the key experimental achievements and conclusions regarding the self-assembly of spherical 
colloids via DNA-mediated interactions I17II23I . 

Conjugation and characterization methods, three representative self-assembly strategies of the self- 
assembly of DNA NPs into super-lattices O, and design principles for plasmonic nano-structures I10II331 - 
|35) have been recently discussed, for example, how DNA has been used to build finite-number assemblies 
(plasmonic molecules), molecularly driven plasmonic switches (36) . regularly spaced NP chains (plas¬ 
monic polymers), and extended two- (2D) and three-dimensional (3D) ordered arrays (plasmonic crys¬ 
tals) (To) [^, as well as DNA patchy particles (33 and oligonucleotide capsules (^, and the formation 
and stability of surface-tethered DNA gold-dendron NPs (43 have been discussed. Some work has dealt 
with the plasmid-templated shape control of condensed DNA block copolymer NPs (41). Of particular in¬ 
terest has also been the integration of self-assembly with top-down lithography (21) for nano-patterning 
of materials onto substrates (42). 

Geerts et al. (23 have reviewed various aspects of DNA-functionalized NPs self-assembly, such as the 
synthesis of DNA NPs nano-complexes, their physical properties and a range of applications for nano- 
and micron-sized colloids either in two- or three-strand systems. At the same time, in two-strand systems. 
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the DNA strands belonging to different NPs hybridize directly to each other, in three-strand systems a 
third DNA linker is required for binding the DNA strands belonging to different NPs. The effect of the 
number of base pair matches, length, DNA grafting density on the colloidal surface, ionic strength of 
the surrounding medium and effects of pH and solvent polarity, etc., are parameters discussed in this 
work. In particular, understanding and controlling DNA-mediated interactions between colloidal parti¬ 
cles is a key to the self-assembly process. A number of factors contributing to these interactions have 
been investigated by experiments, theoretical studies and computer simulations designed to investigate 
interactions and aggregation/melting behavior of DNA-functionalized NPs (^. Moreover, the clustering 
of DNA-functionalized colloids that occurs during the self-assembly shares similarities with percolation 
phenomena {4^. Geerts et al. (20) have also discussed ID, 2D and 3D assemblies of DNA-functionalized 
colloids, where for the 3D assemblies an account of micron- and nano-sized NPs is given. Additionally, a 
discussion on the broad range of applications of such systems is provided focusing on applications such as 
the detection of small molecules with DNA sequences (bio-recognition, e.g., within the cell), or DNA/RNA 
delivery into cells (20) . 

The self-assembly of DNA-functionalized coUoids has been interesting in a number of different 
contexts apart from the self-assembly. For example, recent advances in the synthesis of bio-molecule 
NP/nanorod hybrid systems applications in the generation of 2D and 3D ordered structures in solutions 
and on surfaces have been a large area of research (22). The potential use of bacteria as a novel biotech¬ 
nology to facihtate the production of NPs is also an important application, while emphasis has been also 
put on bio-molecule-NP assemblies for bio-analytical applications and for the fabrication of bio-electronic 
devices (22). Examples of such applications are sensors, bio-molecule-functionalized magnetic particles, 
and bio-molecule-based nano-circuity and viral templates (22). Also, composite assemblies of nucleic 
acids, proteins and NPs have been discussed in the literature (22), as well as the protection and promotion 
of UV radiation-induced liposome leakage via DNA-directed assembly with gold (Au) NPs I44II45I . 

Heuer et al. Q highlighted the importance of Au functionalized conjugates for biomedical applica¬ 
tions (e.g., biological sensing (46)) while suggesting a copper-free click chemistry for the programmed 
ligation of DNA-functionalized Au-NP (47) . The focus has been on the invention and applications of “nano¬ 
flares” 14811501 . which carry reporter duplexes being released after binding. This signal provides the abil¬ 
ity to detect specific endo-cellular targets such as mRNAs, microRNAs, or ions and small molecules in real 
time. Their properties, stability, ceUular uptake and cytotoxicity have been the main focus for such sys¬ 
tems. However, the design of many more highly multiplexed systems with improved targeting properties 
are to be expected in the near future. Sokolova et al. (9J discussed the transfection of cell with DNA NPs 
and methods for gene transfer into living cells, while Rosi et al. )8][TT) focused on the use of these nano¬ 
structures in bio-diagnostics, for example electrical and electrochemical detection, magnetic relaxation 
detection, and nano-wire and nano-tube-based detection. Another application relates to the colorimetric 
recognition of DNA intercalators with unmodified Au-NPs (STI . 

Molecular dynamics simulations and theoretical methods related to the self-assembly of DNA- 
functionalized NPs have been briefly discussed previously (^|27][^. A recent review has summarized 
the main theoretical and simulation models in the context of recent experimental progress, highlighting 
the advantages and disadvantages of these models and key findings (29). Moreover, the crystalline struc¬ 
tures obtained so far by experiment have been summarized (29). Despite the above discussion, a com¬ 
plete presentation of the theoretical and simulation contribution to the field is still lacking. This review 
intends to fill this gap and aims at presenting key experimental findings and an overview of theoreti¬ 
cal and simulation studies that have contributed fundamentally to understanding the self-assembly of 
DNA-functionalized NPs. 


2. Experiments 

2.1. Synthiesis and experimental methodology 

A number of DNA-functionalized NPs-based structures can be formed by tuning independently the 
NPs size, shape, and composition. However, a general strategy for the synthesis of DNA-functionalized 
NPs has been discussed only recently (^, where also NPs with high-density shells of nucleic acids can 
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be built. Specific coating technology on the NP surface has been intplemented in order to realize this 
general synthetic strategy of DNA-based self-assentbly. The creation of heterogeneous NP super-lattices 
is based on carboxylic-based conjugation t54l . These heterogeneous materials have novel optical and 
field-responsive properties, which are achieved by a careful selection of NP characteristics, such as size, 
shape, and the number of DNA strands per particle and DNA flexibility An important advance in 

the synthesis of DNA-coated NPs is the synthesis of solid colloids with surface-mobile DNA linkers Il55l . 
In this case, the DNA linkers are fully mobile on the surface of the NP, and the temperature range for 
equilibrium self-assembly is much broader than in the case of immobile linkers. By increasing the tem¬ 
perature window for melting, a better control of the self-assembly process can be achieved, because fast 
crystallization leading to kinetically trapped nano-structures can possibly be avoided. Moreover, link¬ 
ers accumulate at the NPs below the melting temperature, thus increasing the number of hybridization 
events and the stability of the final assembly structure (^. A recent computer simulation study has 
drawn similar conclusions 1^ . 

Although initially the self-assembly was based on Au NPs, self-assembly based on silver NPs conju¬ 
gates by using triple cyclic disulfide moiefies (53 or metal oxide NPs (^ has been also reported. DNA- 
induced size-selective separation of mixtures of Au-NPs is also possible by using the co-operative 
binding properties of DNA-functionalized NPs. This is based on the fact that the melting temperature 
of the NPs aggregates depends on the size of the NPs (60l. Then, a combination of sedimentation and 
centrifugation techniques can be used to separate larger NPs from smaller ones (^. This method of sep¬ 
aration of NPs according to their size has also been demonstrated in the case of a ternary mixture (5^ . 
Controlled structural evolution of large silver NPs and their DNA-mediated bimetallic reversible assem¬ 
blies (61] has been recently studied by Kim et al. (63, as well as the synthesis and reversible assembly of 
DNA Au-NP cluster conjugates (63. 

It has been recently shown that in binary mixtures of sub-micron colloids, certain crystals can be 
obtained from parent crystals by diffusionless transformations, analogous to the Martensitic transforma¬ 
tion of iron (64l . This also provides a possibility of accessing certain structures by these parent crystals, 
which otherwise would be kinetically unaccessible, though this structure has a lower energy. An exam¬ 
ple of this approach is the formation of a CuAu-I structure from a CsCl crystal structure (64) • This work 
underlines the fact that structural transformation of atomic solids may also occur in DNA-linked NPs. An 
exciting feature of DNA NP assembly structures is that the NP crystals form reversibly during heating and 
cooling cycles. For example, the bcc lattice structure is temperature tunable and non-compact with par¬ 
ticles occupying only a small fraction of the volume of the unit cell (63[66j. The compression properties 
of DNA NPs assemblies have been studied by measuring their response to the applied osmotic pressure 
showing that lattices based on these materials are super-compressible and capable of transforming be¬ 
tween different structures (13. This explains the capability of these systems of transforming between 
different structures. The lattices of nano-particles interconnected with DNA, exhibit an isotropic trans¬ 
formation under compression with a remarkably strong decrease of the lattice constant, up to a factor of 
about 1.8, corresponding to more than 80% of the volume reduction (^. The force field extracted from 
experimental data can be well described by a theoretical model that takes into account the confinement 
of DNA chains in the interstitial regions. It has been shown that compression properties of these systems 
can be tuned via DNA molecular design (63- Moreover, the structural plasticity of bio-molecules and the 
reversibility of their interactions can also be used to make nano-structures that are dynamic, reconfig- 
urable, and responsive. The inter-particle distances in the super-lattices and clusters can be modified, 
while preserving structural integrity by adding molecular stimuli (simple DNA strands) after the self- 
assembly processes have been completed. Both systems were found to switch between two distinct rigid 
states, but a transition to a flexible state configuration within a super-lattice has shown a significant hys¬ 
teresis (68) . 

The crystal growth depends on the length of DNA linker, the solution temperature, and the self¬ 
complementarity of DNA strands and can be compared to the case where non-self-complementary DNA 
linker strands are used (63. A three-step process for crystal growth has been suggested, i.e., an initial 
“random binding”, resulting in disorder Au-NP aggregates, a localized reorganization, and a subsequent 
growth of crystalline domain size, where the resulting crystals are well-ordered at aU subsequent stages 
of growth (63- A stepwise growth process to systematically study and control the evolution of a bcc crys¬ 
talline thin-film of DNA-functionalized NPs on a complementary DNA substrate has been recently pre- 
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sented (TOj. This crystal growth on a substrate has been examined as a function of the temperature, the 
number of layers, as well as the substrate-particle bonding interactions. In this case, the interfacial energy 
between various crystal planes and the substrate were tuned by the DNA interactions controlling the crys¬ 
tal orientation and size in a stepwise fashion using chemically programmable attractive forces. This is a 
unique approach, since the prior studies involving super-lattice assembly typically relied on repulsive in¬ 
teractions between particles to dictate the structure, and those that relied on attractive forces (e.g. ionic 
systems) still maintained repulsive particle-substrate interactions (TOl. Indeed, to date, 17 unique sym¬ 
metries have been realized and over 100 unique crystal structures have been synthesized, aU of which 
conform to a key hypothesis: the obtained assembly structures maximize the total number of hybridized 
DNA interconnects between particles (70). 

The number of possible lattices for synthetically programmable NP super-lattices can be increased 
by using a hollow 3D spacer approach (71). In this approach, some NPs within a binary lattice can be 
replaced with “spacer” entities that are designed in such a way as to replace the behavior of the miss¬ 
ing NPs. For example, the NPs of one type are omitted by the structure without affecting the positions of 
the other NPs. In this way, super-lattices with five distinct symmetries have been created, including the 
one that has never been observed before in any crystalline material (71), while the nano-meter precision 
in lattice parameter tunability has been maintained (71). Although hexagonally close-packed arrays of 
NPs are quite commonly made by different assembly strategies, simple hexagonal lattices of NPs have 
been more difficult to create, because close-packed structures are often energetically more favorable. 
Furthermore, a graphite arrangement of spherical NPs has never been made by any assembly method. 
However, by applying the 3D spacer approach to the AB 2 lattice, both of these lattices are easily synthe¬ 
sized and exhibit a crystalline order with micro-meter-sized domains (71). Hollow triangular gold-silver 
alloy nano-boxes have been also reported (72), a process which is reversible upon heating with sharp 
melting transitions, which has been monitored at wavelengths throughout the visible and near-IR. 

2.2. Phase behavior and influence of various parameters 

An elegant way to regulate the range and magnitude of forces between pairs of different particles is 
based on the control of the specific DNA linkages between colloidal particles self-assembly through the 
combination with polymer brushes (^. Such self-assembly has been found to be a reversible process. 
The key of reversibility is not to allow the particles to be very close, where van der Waals forces become 
important. This is why the presence of adsorbed polymers is required (7^. In fact, DNA-functionalized 
Au NPs in macromolecularly crowded polymer solutions have been recently reported (74). 

The inter-particle distance of fee structures can be tuned by increasing the DNA linker length or the 
total DNA length (75). Structures containing less than 0.52% of inorganic material have been synthesized, 
enabling the tuning of the unit cell lattice parameters within a large range of values (^. Moreover, the 
stability of the structures depends on the sequence of the DNA bases (76l[73. Namely, different DNA se¬ 
quences will guide the self-assembly of these NPs into different structures (78). In binary mixtures, the 
degree of binding between complementary beads depends also on the number of matching pairs and 
the ionic strength of the solution (79). A large variety of colloidal structures, such as chains of alternat¬ 
ing large and small particles for a large range of particle number ratios and volume fractions have been 
investigated (^. Also, the plasmon coupling strength depends strongly on the length of the DNA connect¬ 
ing Au-NPs go). The choice of the DNA linking molecules and the absence or presence of a non-bonding 
single-base flexor can be adjusted, so that Au NPs assemble into micro-meter-sized fee or bcc crystal 
structures. These results demonstrate that a single-component system can be directed to form different 
structures (78). Moreover, particles with a poly-dispersity of 20% do not form weU-defined crystalline 
assemblies. Crystalline structures usually require a poly-dispersity of NPs below 10% (78). 

For 3D hybridization with polyvalent DNA Au NP conjugates (iQ, it has been found that there is a 
minimum number of base pairings necessary to stabilize DNA Au NP aggregates as a function of salt 
concentration for particles between 15 and 150 nm in diameter. Sequences containing a single base pair 
interaction are capable of effecting hybridization between 150 nm DNA Au NPs. While traditional DNA 
hybridization involves two strands interacting in one dimension (ID, Z), hybridization in the context of 
an aggregate of polyvalent DNA Au NP conjugates occurs in 3D making NP assembly possible with three 
or fewer base pairings per DNA strand. These studies enable the comparison of the stability of duplex 
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DNA free in solution and bound to the NP surface. It has been found that 4-8, 6-19, or 8-33 additional 
DNA bases must be added to free duplex DNA to achieve melting temperatures equivalent to hybridized 
systems formed from IS, 60, or ISO nm DNA Au NPs, respectively (16). In addition, the equilibrium bind¬ 
ing constant (K^q) for IS nm DNA Au NPs (three base pairs) is three orders of magnitude higher than the 
Keq of the corresponding NP free system (16). 

Xiong et al. )8T][82) provided the first experimental phase diagram of the NPs assembled by DNA link¬ 
ers in terms of nominal DNA linker/Au-NP ratio r and volume fraction rj, where the boundary between 
disordered and bcc structures has been determined. The formation of a crystalline bcc phase was ob¬ 
served for a hmited range of linker lengths, while the number of linkers per particle controlled the onset 
of system crystallization. The effect of linkage defects on the crystalline structure was also examined. 
Note that these results were corroborated later by computer simulations (83l . 

Although a broad range of different non-compact structures has been created over the last years (84), 
a single-component NP assembly of a diamond structure has not been reported. By a delicate tuning of 
the particle interactions, binary diamond-like crystalline structures have been fabricated ^4) though. For 
example, NaTl-type non-compact structures have been reported, using surface-modified Qp phage capsid 
particles and Au-NPs of similar radii. This structure has an inorganic component and an organic virus-like 
protein NP interpenetrating diamond lattices. This is an example where particles of different chemical 
nature can be combined after a proper surface processing, without phase separation taking place. How¬ 
ever, protein-based structures had been discussed already before (85), in addition to the inorganic NPs 
(86). NP-DNA conjugates bearing a specific number of short DNA strands by enzymatic manipulation of 
NP-bound DNA have been also reported (87), as well as the self-assembly of icosahedral virus particles 
organized by attached oligo-nucleotides (88), and binary heterogeneous super-lattices assembled from 
quantum dots and Au-NPs with DNA (83, or assembling colloidal clusters using crystalline templates and 
reprogrammable DNA interactions 190) . 


2.3. Micron-sized nano-particles 

Micron-sized latex particles with single-stranded DNA grafted to the surface have been used as a 
model system to study DNA mediated interactions (91). The presence of DNA strands allows for the tun¬ 
ing of the interactions between NPs by combining hybridizing linker DNA with non-hybridizing neutral 
DNA, which leads to a gradual change of the assembly rate. The effect of linker/neutral DNA ratios on 
particle assembly kinetics and aggregate morphology has been experimentally investigated for a range 
of ionic strengths. The conditions for controlling various assembly morphologies have been identified, 
and the involved attractive and repulsive interactions have been described and explained for the pro¬ 
posed approach, also accounting for further theoretical considerations (91). 

Kim et al. (92]|93) have reported on the first DNA-mediated micrometer sized colloidal crystals by us¬ 
ing polystyrene particles, discussing preparation strategies of these coUoids in order to obtain such crys¬ 
tals. Thermodynamic properties of the self-assembly of these crystals are also predicted through a simple 
thermodynamic model. Possible ways of obtaining more ordered structures by studying the growth ki¬ 
netics and phase behavior of the self-assembly of DNA-functionalized colloids have been also considered 
1^ . Various aspects of the selective, controllable, and reversible aggregation of polystyrene Latex micro¬ 
spheres via DNA hybridization have been discussed (94)- These systems were found to have interesting 
kinetics and non-exponential binding behavior (95) . 

The first direct measurements of DNA-induced interactions between a pair of micron-sized polymer 
colloids have been conducted by Biancaniello et al. (^ by using an optical tweezer method. These re¬ 
sults can be modelled using standard statistical mechanics methods. It has been shown that micron-sized 
spherical colloids have binding dynamics of a power-law in annealing and crystallization processes. Ad¬ 
ditionally, the phase behavior of DNA-mediated micro-sphere suspensions was determined as a function 
of salt and oligonucleotide concentration (93. In this case, fluid phases, aggregate phases, or mixed fluid 
phases with aggregates have been observed, while increasing salt or hybridizing events favors the forma¬ 
tion of aggregates. Forces of interaction between DNA-grafted colloids via optical tweezer measurements 
have been also reported elsewhere (98), in quantitative agreement with a proposed theoretical model 
(99). The specific inferaction of DNA-functionalized polymer colloids by using DNA linkers in solution has 
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been also reported llOOl . Based on such measurements, a Monte Carlo (MC) scheme to describe colloidal 
crystallization has been implemented I101II102I . 

Switchable self-protected attractions in DNA-functionalized spherical colloids are also very interest¬ 
ing 11031 . For example, single-stranded DNA can fold creating secondary structures, such as loops and 
hairpins I104itl071 . which influence the self-assembly. For micron-sized particles, these structures pro¬ 
vide a control over inter-particle binding strength and association kinetics. By heating and cooling the 
system, one can mediate these effects. Such topological interactions may lead to new materials and phe¬ 
nomena, like particles strung on necklaces, confined motions on designed contours and surfaces, or even 
colloidal arrangement in the shape of “Olympic gels” 11041 . Also, polyvalent nucleic acid nano-structures 
(PNANs), which are composed of only cross-linked and oriented nucleic acids, are of great interest. Apart 
from the self-assembly properties, these particles are capable of effecting high cellular uptake and gene 
regulation without the need of a cationic polymer co-carrier llOSl . 

2.4. Anisotropic DNA-functionalized nanoparticles 

Some structures cannot be fabricated by using only spherical colloids in all dimensionalities. Hence, 
directional bonding interactions which can be tuned by the different geometry of the NP have been used 
for a better control over the directed self-assembly of such DNA NPs. One advantage of using particle 
anisotropy is the fact that a greater surface contact favors a higher number of hybridizations. More¬ 
over, the effective local concentration of terminal DNA nucleotides that mediates hybridization increases, 
while conformational stresses imposed on nano-particle-bound ligands participating in interactions be¬ 
tween curved surfaces might be smaller iFTl 11091 . DNA NPs with core particles of different shape have 
been synthesized in a way to study the effect of the core-core particle interactions, affecting super-lattice 
dimensionality, crystallographic symmetry and phase behavior I109II110I . 

An example of ID structures is the self-assembly of nano-rods Hill , where again the most stable 
structure is the one that maximizes hybridization events between DNA sticky ends (2^. Moreover, for 
octahedron NPs, one might obtain fee or bcc structures depending on the length of the DNA strands. 
However, previous work 11091 IllOl has obtained self-assembly aggregates rather than real long-range 
crystal structures indicating the difficulty in obtaining long-range crystalline structures. Recently, the use 
of DNA-drlven assembled phospholipid nano-discs as a scaffold for Au-NP patterning has been discussed 
11121 . A facile and efficient preparation of anisotropic DNA-functionalized Au NPs and their assembly has 
been reported recently 11131 . This one-step solution-based method has been used to prepare anisotropic 
DNA-functionalized Au-NPs with 96% yield and with high DNA density. Well-defined nano-assemblies us¬ 
ing Au NPs having a specific number of DNA strands have been also discussed 11141 . Periodic square-like 
Au NP arrays templated by self-assembled 2D DNA nano-grids on a surface have been also reported in the 
literature 11151 . as well as multiple-nano-component arrays by 2D DNA scaffolding 11161 . Triangular gold 
nano-prisms 11171 in the presence of attractive depletion forces and repulsive electrostatic forces assem¬ 
ble Into equilibrium ID lamellar crystals In solution with inter-particle spacings greater than four times 
the thickness of the nano-prisms IllSl . Experimental and theoretical studies reveal that the anomalously 
large spacings of the lamellar morphologies are due to the balance between depletion and electrostatic 
interactions 11181 . 

DNA-directed assembly of asymmetric nano-clusters using Janus NPs has been also investigated 
I119II120I . Curvature effects in DNA Au-NP conjugates are notably important, due to the asymmetry of the 
NP I121II122I . Particle curvature also plays an important role In controlling the DNA surface density 11221 . 
Based on such an approach, there is a method to predict DNA packing on non-spherical particles. An ex¬ 
ample of usage has been provided for Au nano-rods 11221 . Linear micro-structures in DNA nano-rod self- 
assembly have been reported recently 11231 . Fibers, tubules, and ribbons are typical ID nano-structures 
that require anisotropic interactions that also emerge from the geometry of the NP, which affects the 
underlying inter-particle Interactions during the structure formation. For systems of DNA-functionalized 
nano-rods, the resulting structures are ID ladder-like meso-scale ribbons with a side-by-side rod arrange¬ 
ment. Moreover, the DNA can bind reversibly, facilitating in this way the formation of hierarchical as¬ 
semblies with time. Linear structures of alternating rods and spheres have been also obtained, implying 
the generic applicability of the mechanism for nano-scale objects interacting via flexible multiple linkers 
112 ^ . 
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2.5. Self-assembly in one- and two-dimensional geometries 

DNA-based self-assembly of NPs has been discussed in the context of planar or other geometries, 
as, for example, in the case of DNA-mediated 2D colloidal crystallization above different attractive sur¬ 
faces 11241 . Cheng et al. 11251 reported the first DNA-based strategy for mono-layered free-standing NP 
super-lattices or very thin films. In this case, discrete free-standing super-lattice sheets can be obtained 
without the requirement of Watson Crick’s specific binding in which inter-particle spacing occurs within 
the structrure, leading to the formation of plasmonic and mechanical properties that can be rationally 
tuned by adjusting the DNA length 11261 . Moreover, thin standing films based on DNA-coated colloids with 
double-stranded DNA have been reported 11271 . Such colloids can self-assemble into unique crystalline 
mono-layers, suspended at a distance of several colloidal diameters above a weakly adsorbing substrate, 
without a dependence on DNA hybridization. The advantage of such systems is that they can be prepared 
in one location and then be used in a different application somewhere else. This may be a way for the 
assembly of multi-component, layered colloidal crystals 11271 . Moreover, responsive multi-domain free¬ 
standing films of spherical Au-NPs assembled by DNA-directed layer-by-layer approach have been also 
studied recently 11281 . Large-area spatially ordered arrays of Au-NP crystals directed by lithographically 
confined DNA origami have been also synfhesized 11291 and 5 nm size Au-NPs arranged info long-range 
ordered 2D arrays onfo lithographically patterned substrates Il29l . Additionally, the first DNA-mediated 
formation of supra-molecular mono- and multi-layered NP structures has been presented 11301 . as well 
as DNA and DNAzyme-mediated 2D colloidal assembly 11311 . 


3. Computer simulations 

Various simulation models for DNA chains have been contrived over the last years 1132111431 . As it 
usually happens with these models, they aim at reproducing only the properties of interest of the system. 
Due to the large system sizes and time scales required in the study of DNA or DNA-coated NPs, mainly 
coarse-grained models are capable of describing relevant phenomena or bespoke theoretical models, 
such as the self-assembly of DNA-functionahzed NPs. 

3.1. Coarse-grained models 

Starr and Sciortino I142I have proposed a simple model for the study of the self-assembly of DNA- 
functionalized NPs (figure [^. This model has been used for short DNA strands (six to twelve bases) and 
specific DNA sequences. The effect of several parameters has been explored, for example, the number of 
bases per strand, the number of linking bases and the number of spacer bases on the stability of crystal 
states, as well as the effect of using a single linking NP t)q)e versus a binary linking system. Additionally, 
the free energy, entropy, and melting point for bcc and fee lattices have been explicitly calculated. While 
binary systems preferably form bcc lattices, melts of single NP type with chains of self-complementary 
basis form the most stable fee crystals. A way to maximize a crystal stability based on the heat of fusion 
between crystal and amorphous phases growth has been discussed I1441 . The stability of bcc and fee 
crystals formed by DNA-linked NPs has been further investigated I145I as a function of the number of 
attached strands to the NP surface, the size of the NP core, and the rigidity of the strand attachment. It has 
been found for this model that relaxing the condition of kineticaUy immobile strands on the surface of the 
particles, a slight increase or decrease in the melting temperature occurs. Larger changes to the melting 
temperature result from increasing the number of strands, which increases the melting temperature, 
or by increasing the core NP diameter, which decreases the melting temperature ll43l . A quantitative 
description of the melting temperature has been provided by the model of Starr and Sciortino 1142111441 . 
This model has been fully investigated and used in the study of self-assembly process I146II147I . In these 
studies, a phase diagram for the NPs as a function of temperature and NPs density has been provided. 

Furthermore, the effective potential between DNA-functionalized NPs can be obtained enabling the 
comparison of the phase diagram of effective particles with the phase diagram of the NPs based on the 
fuU coarse-grained description. Various grafting geometries have been considered in studying the self- 
assembly, such as triangular planar, tetrahedral, square-based pyramid, triangle-based bipyramid, and 
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Figure 1. (Color online) Schematic representation of the model of Starr and Sciortino 11421 . This model 
considers only a specific DNA sequence (specified in panel a) which leads to the hybridization shown 
in panel (b). This model can assume different DNA lengths, NPs sizes and grafting geometries (c)-(e). 
The effective potential between two NPs modelled in this way depends on the angles 6i and 02. and the 
distance r between the NPs. Further details on how to calculate the effective potential can be found in 
11431 . Figure adapted from 11431 after permission. 


octahedral geometries 11471 . These NPs have short DNA strands with imposed rigidity and complemen¬ 
tary DNA sequences resembling interactions between patchy colloids. For all these geometries, but the 
octahedral one, there has been found a phase behavior analogous to gas-liquid transition. For the octahe¬ 
dral geometry, various crystal structures of similar symmetry occur, depending on the chemical potential 
(density) of the system. Then, interpenetrating networks of DNA-functionalized NPs occur 11471 . In the 
following study Il48l . a two-step crystallization of DNA-functionalized NPs has been discussed, where 
the narrow temperature range that the self-assembly takes place is described through a two-step process, 
mediated by a highly connected amorphous intermediate. Moreover, at lower temperatures, the system 
is kinetically trapped hindering crystallization 11481 . By using a standard cluster analysis and relevant 
order parameters based on spherical harmonics, the structures and their degree of crystaUinity were re¬ 
solved. The relaxation and crystallization times of such systems have been also studied Il48l . Different 
amorphous networked phases induced by multiple liquid-liquid critical points providing a phase diagram 
of the temperature versus density have been identified 11491 . based on effective NPs. This effective po¬ 
tential has been carefully revised in the following work 11431 . Different structure factors can describe the 
interpenetrating networks as the density of the NPs increases 11491 . 

The equilibrium clustering dynamics and the self-assembly kinetics can be described on the basis of 
the Flory-Stockmayer theory 1150111521 finding good agreement among theory and results obtained from 
the coarse-grained model. Considering only two- and three-fold NPs in order to facilitate the branching 
between two-fold DNA NPs, percolation and phase separation at low density have been studied 11521 . 
Furthermore, the kinetics of the self-assembly are investigated in terms of the bond formation rate as a 
function of temperature and time 11521 . It has been found that these processes follow an Arrhenius behav¬ 
ior, where the self-assembly is described by two “reaction” constants. The model of Starr and Sciortino 
11421 has been also used to study the effects of bidispersity in DNA strand length 11531 and effects of 
DNA strand sequence and composition I153II154I . Moreover, NP dimers linked by DNA can be used as 
the building blocks of a multi-scale, hierarchical assembly, especially when the size of the NPs is similar 
11551 . The structure of NP dimers through a combination of scattering experiments and molecular sim- 
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ulations has been recently investigated 11551 . It has been found that the inter-particle separation within 
the dimer is mainly controlled by the number of hybridized DNA strands, their length and the curva¬ 
ture of the particle, as well as the excluded volume effects. Without any free parameters, the increase of 
dimer separation with increasing temperature can be interpreted as a result of the change in the number 
of connecting DNA 11551 . Recently, it has been shown that tetravalent DNA nano-stars never crystallize, 
but they form a fuUy bonded equilibrium gel 11561 . The phase diagram of these organic complexes has 
been also discussed 11571 . 



Figure 2. (Color online) Schematic representation of the model by Li et al. 11611 . as adapted from the 
original model of Knorowski et al. (83) . In panel (a) we present a typical DNA chain. Blue beads indi¬ 
cate flexible ssDNA, red beads dsDNA, while green, black, and orange beads represent DNA bases. The 
magenta flankening beads have been originally used in (83) to prevent double bonding, in this way over¬ 
coming a drawback of the model of Starr and Sciortino 11421 . Li et al. have used spherical particles for 
their simulations, but different geometries can be used as those in panels (b), (c), and (d). In panel (e), we 
present an initial configuration for MD or MC simulation of DNA-functionalized octahedra based on the 
model of Li et al. 11611 . 

A model that is closer to the experimental expectations has been proposed by Knorowski et al. 
I83II158|[I5^ and has been amended by Li et al. 11601 to match the experimental energy and length scales 
(figurej^. Apart from a number of different structures and the very good prediction of the experimental 
phase diagram as a function of the volume fraction, temperature, and grafting density, this model is capa¬ 
ble of describing the dynamics and thermodynamics of the self-assembly process. The D-bcc crystalline 
phase also exists in their phase diagram. In this case, two different types of DNA-functionalized NPs oc¬ 
cupy the positions of a bcc lattice, but the type of particles does not have the same pattern across the crys¬ 
talline structure. This effect has been previously overlooked. Moreover, a CsCl-bcc structure is a part of 
the phase diagram, while the exact boundary with disordered structures has been determined, providing 
the universal dynamics of the crystallization process (83) . This model has in many different aspects over¬ 
come the problems of the modelling approach of 11421 . such as the possibility of a strand to bind to more 
than one complementary DNA strands. Recently, this model has been used to study in more detail the crys¬ 
tallization dynamics, such as gelation, nucleation, and topological defects 11591 . The formation of crystals 
involves various dynamic stages of the system, starting from a liquid mixture to the formation of a gel. 
Then, the formation of a crystal containing defects takes place, and finally an equilibrium structure forms 
11591 . For these stages, the dynamics of the process have been discussed in detail Il59l . Li et al. 11601 have 
actually used this model (83) to study the phase diagrams as a function of size ratio and DNA coverage 
ratio as is done in experiment (^. The authors confirmed fully the experimental observations I23II31I . 
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In a follow-up study 11611 . it was found that the therntally active hybridization drives the crystalliza¬ 
tion of DNA-functionalized NPs. By using Molecular Dynamics (MD) simulations, the dynamic aspects of 
the assembly process were analyzed and the key ingredients to a successful assembly of NP super-lattices 
have been identified. The conclusions of the simulation model are scalable, being able to reproduce the 
experimental observations, which have proposed various structures 11611 . Moreover, the optimum per¬ 
centage of DNA hybridization events has been determined for different binary systems, being able to 
propose suitable linker sequences for future nano-material design 11611 . Recently, this model, in combi¬ 
nation with experiments, has been used to investigate the DNA-mediated NP crystalhzation into Wulff 
polyhedra, where the use of DNA-functionalized NPs offer a higher control over the growth process of 
the crystals 11621 . 

The pair interactions between DNA-functionalized colloids have been the main topic in the field 11371 . 
and among other methods, MC simulations considering NPs with short and stiff DNA chains have been 
employed 11631 . Such interactions have been initially measured for planar surfaces corresponding to 
large NPs. Based on simulation results, analytical expressions for the interactions can be obtained, claim¬ 
ing that this is a valid description for a range of NPs from the nano- to the micro-meter scales, as well as 
for a large range of grafting densities. Moreover, important contributions to the repulsive and attractive 
terms of the free energy are coming from purely entropic effects of the discrete tethered sticky ends. 
Such entropic contributions compare well with the hybridization free energy of a free pair of sticky ends 
in solutions. Stable gas-liquid separation only occurs for NPs with radii of a few tens of nanometers, sug¬ 
gesting different crystallization routes for nano- and micron-sized NPs Il63l . Statistical anisotropy on the 
distribution of sticky ends on the NPs’ surface leads to a large variation of the binding strength. This 
may weaken the reliability of tests based on the detection of DNA targets in diagnostics. These results 
provide a general background for the systems with tethered binding groups 11641 . as for example in the 
areas of supra-molecular chemistry or ligand/receptor mediated bio-recognition (l63l . Based on this MC 
coarse-grained model, a new strategy to improve the self-assembly properties of DNA-functionalized col¬ 
loids has been proposed 11651 . The narrow crystallization temperature window is due to the sensitivity 
of DNA-mediated interactions to temperature and other physical control parameters. By exploiting the 
competition between DNA linkages with different nucleotide sequences, this temperature window can 
become larger. Thus, the temperature dependence of the self-assembly process has been decreased offer¬ 
ing a control over a larger temperature range. Such an approach is also applicable in a realistic situation, 
where the control is due to a combinatorial entropy gain, which stems from the large number of bind¬ 
ing partners for each DNA strand. Experimental work has suggested particular DNA sequences that are 
capable of reproducing those effects 11651 . 

The narrowness of the temperature has been also addressed in another study by Uberti et al. |166| . 
since this narrow temperature window from dilute systems to aggregates reduces considerably the ca¬ 
pability of obtaining defectless crystal structures, thus restricting the number of possible crystal struc¬ 
tures. A design of DNA-coated colloids which crystallize on cooling and melt on further cooling has been 
proposed 11661 . Although this is a mere theoretical model-assumption, it has been argued that realistic 
situations that can follow such a phase behavior do exist. While the results presented refer to binary sys¬ 
tems, these concepts also apply in multi-component systems, opening the way towards the design of truly 
complex self-assembling structures based on spherical colloids Il66l . The imphcations of these results 
have been commented by Gang 11671 . Latest techniques exploit the narrow temperature window though, 
in order to construct amorphous materials with specific morphology and local separation between mul¬ 
tiple components. In this regard, strategies to direct the gelation of two-component colloidal mixtures 
by sequentially activating selective interactions have been proposed 11681 . Changes in the structure can 
be investigated via MD simulations and experimental methods. This approach in studying the multi-step 
kinetics of the self-assembly process can be extended to the assemblies of multi-component meso-porous 
materials with possible applications in hybrid photovoltaics, photonics, and drug delivery 11681 . 

MC simulations combined with histogram re-weighting techniques have been applied in order to con¬ 
struct the phase diagram of DNA-functionalized NPs based on a coarse-grained model 11691 . This model 
can predict the phase separation between a dilute vapor-like phase and a dense phase where NPs are a 
part of a liquid-like network based on DNA links. In this case, there is a non-monotonous dependence of 
the coexistence pressure on the temperature, namely the reduced hybridization free energy 11691 . This 
behavior may be attributed to the crossover between two distinct regimes. Namely, the hybridization- 
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free-energy driven regime and an entropy-driven regime. In the former case, there is a “normal” vapor- 
liquid equUibrium, where the system gains free energy but loses entropy. In the entropy-driven regime, 
an increase in entropy occurs due to the rearrangement of sticky-end bonds in the liquid phase. The main 
conclusion of this work is that the system undergoes a phase separation if the number of DNA-chains per 
NP is larger than two. As the number of chains increase, the coexistence region increases considerably 
11691 . In another MC study 11701 . the phase behavior of coUoids coated with long, flexible DNA chains, 
which can also be achieved experimentaUy 11711 . has been investigated. In this case, when the number 
of DNA strands per NP decreases below a critical value, the triple point disappears. Also, the condensed 
phase coexisting with the vapor remains a liquid. This result may explain the reason why NPs with small 
DNA coverage cannot crystallize in dilute solutions highlighting in this way the discrete nature of DNA 
binding 11701 . 

In another coarse-grained approach realized in a MC scheme, a simplified model was proposed, based 
on experimental data, which can describe the phase behavior for DNA-coated colloids 11721 . These results 
have indicated that the assumption of pair-wise additivity of DNA-mediated interactions lead to substan¬ 
tial errors in the estimation of the free energy of the crystal phase. While these simulations correctly 
predict the experimentaUy obtained crystal structures, despite uncertainties in the estimates of contour 
and persistence lengths of ssDNA chains, they also provide a good estimate of the melting temperature 
11721 . Hence, further, more accurate experimental results are required in order to construct a more ac¬ 
curate coarse-grained approach 11721 . 

The effect of linear shear force on the self-assembly of DNA-functionalized colloids in two dimensions 
has been also investigated 11731 . A better sampling of the phase space is possible in binary systems. From 
the constructed phase diagram, a linear dependence of the minimum DNA force necessary for the line 
formation on the dimensionless shear rate was found. By employing a force analysis on the obtained 
structures, a linear orientation perpendicular to the axis of the elongations component of the shear was 
found, which is the orientation that enables the DNA attraction to resist to shear 11731 . The effect of 
internal viscoelastic modes on the Brownian motion of a 1-DNA coated colloid is yet another topic of 
recent research Il74l . 

Spherical NPs of 1 mm colloids bridged by DNA with 32mm contour length have been studied in a 
mixture of two types of particles with grafted double-stranded 1-DNA with short complementary single 
stranded overhangs as free binding ends. Confocal microscopy identified the formation of stable size- 
limited clusters in which the two types of particles are in contact. The simulations suggest that entropic 
exclusion of the bridging DNA from space between the nearby particle surfaces hinders the formation of 
crystals in milimeter-sized NPs 11751 . 

3.2. Further models, theoretical, and numerical approaches 

Specific attraction between NPs can be induced by DNA bridging creating a large number of peri¬ 
odic structures. Although the thermodynamics of DNA-functionalized NPs in solution is well-understood, 
perfectly suitable coarse-grained models are still challenging. This difficulty has hindered the design of 
complex temperature, sequence and time-dependent interactions required to describe, for example, ma¬ 
terials with highly complex or multi-component micro-structures or the abihty to reconfigure or self- 
replicate. A recent study ll76l has provided high resolution measurements of DNA-mediated interactions 
in polystyrene micro-sphere-sized spherical colloids closely related to the experimental conditions. A 
numerically tractable model that quantitatively describes the separation dependence and temperature- 
dependent strength of DNA-mediated interactions can be used, without empirical corrections 11771 . This 
model has been also used to successfully describe the case of grafted DNA brushes with self-interactions 
that compete with the inter-particle formation. This simulation design approach can propose experimen¬ 
tally realizable nano-structures 11771 . despite some criticism of this approach stating that the agreement 
with experiment maybe fortuitous 11781 . 

Similarly, based on standard thermodynamic arguments, a simple quantitative model for the re¬ 
versible association of DNA-coated micron-sized colloids has been built 11791 . These results were able 
to reproduce very well the experimental data regarding the dependence on temperature and sticky end 
coverage indicating clearly the short range of association-dissociation transition of DNA-coated colloids 
I176II179II180I . It has been shown that this transition can be described by a simple model for different 
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types of DNA and different DNA coverage 11791 . This model takes into account the feature of the grafted 
DNA underlining the significance of the entropy cost due to DNA confinement onto the surface of the NPs. 
Here, an extensive description of the experiments carried out to validate this simple theoretical model is 
provided ll79l . Then, a step-by-step statistical mechanics based model is developed, which is proven to 
describe with excellent agreement the temperature and width of the transition, which are both essential 
in the self-assembly of DNA-functionalized colloids 11761 . 

A simple model for the melting and optical properties of DNA Au-NPs aggregates has been also dis¬ 
cussed II 8 II . At the melting temperature, the optical properties of the aggregates change considerably, 
a transition that is narrower for larger NPs, as it has been also discussed elsewhere 11431 . The extinc¬ 
tion coefficient has been computed using the discrete dipole approximation, and an increasing number 
of bonds breaks as the temperature increases. Moreover, the melting temperature approximately corre¬ 
sponds to the bond percolation threshold 11821 . NPs with radius of 50 nm have been considered, which 
compare well with experimentally used NPs. In this MC scheme, there is a reaction-limited cluster-cluster 
aggregation followed by dehybridization of the DNA links. These links form with a probability that de¬ 
pends on the temperature and the particle radius. The final structure is also a function of the number of 
NPs and the relaxation time. At low temperatures, one obtains clusters which after a long relaxation time 
transform into a compact non-fractal cluster. The optical properties I183II184I of these structures can be 
estimated by using the discrete dipole approximation. Despite this structural transformation, the melting 
transition becomes apparent at the extinction coefficient at 520 nm wavelength remaining sharp, while 
the melting temperature increases with increasing NP radius as found in the previous model II 8 II . How¬ 
ever, the restructuring of NPs increases the link fraction at melting to a value well above the percolation 
threshold. The extinction crossover compares well with experiments on Au/DNA nano-composites, where 
the incomplete relaxation is a robust effect. These results share similarities and may relate to a possible 
sol-gel transition in such aggregates II 8 II . Also, photo-switchable DNA-functionalized Au NPs, where hy¬ 
bridization is affected by the “photon dose”, can be used for tuning the DNA-mediated interactions ll85l . 
as well as high-energy Al/CuO nano-composites II 86 I . The optical and topological characterization of Au 
NP dimers linked by a single DNA double strand has been also discussed Il87l . Furthermore, connected 
magnetic micro-particles with DNA G-quadruplexes have been recently reported in the literature 11881 . 

A general theory for predicting the interaction potentials between DNA-coated coUoids has been re¬ 
cently presented I189II190I . This theory incorporates the configurational and combinatorial entropic fac¬ 
tors that play a key role in valence-limited interactions. By enforcing self-consistency, the modelling can 
achieve quantitative agreement to detailed MC calculations. With suitable approximations and in par¬ 
ticular geometries, this theory reduces to the previous successful treatments, which are now united in a 
common and extensible framework 11891 . 

A theory based on an attraction of unlike particles combined with a type-independent soft-core re¬ 
pulsion has been realized for a binary system of DNA-functionalized colloids I19ll . In an experimen¬ 
tal system, the repulsion and attraction may be induced by a DNA solution. Based on a theoretical 
treatment ll91l . the system exhibits diverse and unusual morphologies, such as the diamond lattice and 
the membrane phase with an in-line square order. This is the first time, albeit theoretically, that the di¬ 
amond structure has been predicted Il9ll . Moreover, the inverse approach has been discussed, when a 
particular scheme is proposed for building an arbitrary desired nano-structure. The conditions for ro¬ 
bust self-assembly of the target structures have been identified, which includes the minimum number of 
“colors” needed to encode the inter-particle bonds through the hybridization of DNA pairs. Furthermore, 
it was shown that a floppy network with thermal fluctuations can possess an entropic rigidity and can be 
described as a traditional elastic solid. The onset of the entropic rigidity determines the minimal number 
of bond types per NP needed to encode the desired structure. This theoretical model takes into account 
thermodynamic considerations providing additional conditions for the occurrence of the equilibrium 
structures 11921 . 

Halverson and Tkachenko 11931 have studied the formation of complex mesoscopic structures with 
programmable local morphology and complex overall shape NPs linked with DNA. The hybridization of 
bonds is reahzed through a harmonic potential and a Poissonian statistical distribution. These results 
compare well with experimental results referring to the melting curve of dimers of DNA-functionalized 
NPs, and ordered structures in the range 5-50 nm were obtained. The key elements of this design strategy 
are the symmetry based on the neighbors being of different type and different sequences of bonds, the 
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controlled nucleation, and the cooperative binding. Based on this strategy, a number of complete nano¬ 
structures can be assembled, such as the empire state building 11931 . The size limits of self-assembled 
colloid structures made using specific short-range interactions have been also recently discussed Il94l . It 
has been found that structures with highly variable shapes made out of dozens of particles can form with 
high yield, as long as each particle in the structure binds only to the particles in their local environment 
11941 . Such an approach allows complete enumeration of the energy landscape and gives bounds on how 
large a colloidal structure can assemble with high yield. For large structures, the yield can be significant, 
even with hundreds of particles ll94l . 

A theoretical discussion of a self-assembly scheme which enables the use of DNA to uniquely en¬ 
code the composition and structure of micro- and nano-particle clusters has been presented in a series 
of papers by Licata and Tkachenko I1951I198I . Several important aspects of possible experimental imple¬ 
mentation of the proposed scheme have been discussed. For example, the competition between different 
types of clusters in a solution, possible jamming in an unwanted configuration, and the degeneracy due 
to symmetry with respect to particle permutations Il95l . This methodology has been enabled from the 
high selectivity of DNA-mediated interactions between NPs. Based on a simple model and a more realistic 
description of the DNA-colloidal system, it is possible to suppress disordered glassy phases of the systems 
in order to obtain self-assembly into perfect crystal structures. This requires a combination of stretch¬ 
able inter-particle linkers, which in a realistic situation may simply correspond to sufficiently long DNA 
strands, and, also, a soft repulsive potential. This approach has been implemented on a large number 
of cluster types and is experimentally realizable 11961 . The self-assembly for the case of DNA-caged NPs 
has been also considered in this work 11981 . The dynamics of models that contain such key-lock binding 
interactions, like the Watson-Crick interactions in the DNA including systems where antibodies or cell 
membrane proteins are grafted on NPs can be described theoretically |197| . Depending on the grafting 
density of the binding groups, two regimes have been predicted. In the localized regime, at low coverage, 
the system exhibits an exponential distribution of particle-departure times. At higher coverage, there is 
an interplay between departure dynamics and particle diffusion. This interplay leads to the phenomena 
similar to aging in glassy systems, corresponding to a sharp increase of the departure times. The diffu¬ 
sive effect is similar to dispersive transport in disordered semiconductors. Depending on the interaction 
parameters, the diffusion behavior ranges from standard to anomalous, sub-diffusive behavior 11971 . 

Hydrodynamics may play an important role in the transformation of colloid super-lattice crystallites, 
especially when NPs are micron-sized Il99l . For example, the spontaneous transformation through an¬ 
nealing from a bcc to a fee has been experimentally observed By using computer simulations 

and vibrational mode analysis, one can discern between different structures that have the same energy, 
but hydrodynamic correlations differ. It is the effect of hydrodynamics that induces anisotropic diffusion 
resulting in the formation of different “child” structures during the annealing processes. Therefore, struc¬ 
tures based on DNA-functionalized NPs are not purely determined by thermodynamics, but also require 
the consideration of the process and properties which may be the result of hydrodynamic effects 11991 . 
These phenomena play an important role given the narrow temperature window of the melting curve. 

DNA-functionalized NPs with mobile linkers have been recently studied (56j, inspired by experimen¬ 
tal evidence 155]. This approach combines in a sense the valency of the otherwise isotropic colloids, due 
to the mobility of DNA linkers. By combining theoretical and simulation methods, this valency of colloids 
can be tuned by the nonspecific repulsion between particles. The simulations show that the resulting 
effective interactions lead to low-valency colloids in peculiar “open” structures, considerably different 
from those observed in the case of NPs with immobile DNA linkers (^. 

3.3. Lattice models 

The mapping of microscopic DNA relevant sequences onto the macroscopic phase behavior has been 
in the focus of theoretical and simulation work 1200112021 . Such an approach should enhance the effi¬ 
ciency of methods detecting mutations in specific DNA sequences 12001 . By using a mean-field model, the 
surface and bulk dissolution properties of DNA-functionalized NPs assemblies can be analyzed I201II203I . 
A model for multi-component systems able to form low-symmetry ordered phases has been also pro¬ 
posed despite the spherical symmetry of the interactions 12021 . By combining MC simulations with mean- 
field theory, the thermodynamic, structural and kinetic aspects of stripe phases for simple and multi- 
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component lattice model can be studied 12021 . This lattice model represents a mixture of spherically sym¬ 
metric DNA-functionalized spherical NPs with several species of DNA linkers. However, non-spherically 
ordered structures based on spherical NPs can also be constructed, where DNA linkers of different lengths 
can be used to introduce bonds of different spacing between NPs. The particular example of a stripe struc¬ 
ture has been considered 12021 . In all cases though, the system was forced to have exactly the same ground 
state and the highest degree of connectivity of particles. Hence, as in experiment (23), the most stable 
structure is the one that maximizes the number of hybridization events. However, when the strength of 
competitive binding is weak, a four-component system can exhibit two consecutive ordering transitions, 
which are due to the re-ordering of superstructures within the system. Based on those observations, an 
optimal design strategy has been proposed for super-structures in multi-component systems 12021 . 

Phase transitions in DNA-linked assembhes can be described by a decorated-lattice model 12041 . In 
this study, two types of DNA-linked colloidal mixtures have been considered, namely, systems with iden¬ 
tical NPs functionalized with two different DNA strands and mixtures of two types of particles with each 
one being functionalized with different DNA strands. Establishing the similarity with Ising models with 
temperature and activity dependent effective interactions, various properties for these mixtures can be 
derived, such as the evolution of the dissolution profiles as a function of temperature and number of 
grafted DNA strands, in agreement with experimental systems. The increase of the dissolution temper¬ 
ature with the grafting density is valid only below a certain threshold. For high grafting densities, the 
dissolution temperature becomes a decreasing function of the DNA coverage. In systems with two differ¬ 
ent types of particles, the phase separation involves meta-stable disordered aggregates when compared 
to a phase transition of solvent-rich and ordered phase. The ordering of the coUoidal network enhances 
the stability of the DNA-linked assembly resulting in the dissolution of the aggregates at higher temper¬ 
atures. These results explain the contrasting evolution of the dissolution temperatures with increasing 
probe size in both types of mixtures 1204) . 

A 2D lattice model for particles of many different types with short-range isotropic interactions which 
are pairwise specific has been also introduced 12051 . The main question addressed is whether the ground 
state corresponding to a crystal structure is unique. Interestingly, it has been shown that it is sufficient to 
uniquely determine the ground state for a large class of crystal structures via MC simulations 12051 . 

3.4. Genetic algorithms 

An interesting approach based on genetic algorithms has been recently implemented to design DNA- 
grafted NPs that self-assembly into desired crystalline structures 12061 . Often, the self-assembly first com¬ 
pletes, and then the obtained self-assembly crystal structures are determined. By using genetic algorithms 
based on an experimentally validated evolutionary optimization methodology, it is possible to predict, 
not only original phase diagrams detailing the regions of known crystals, but also several new structures 
(for successful applications of evolutionary algorithms in soft matter systems with directional bonding 
see 1207112101 ). The agreement of these results with the existing experimentally predicted structures val¬ 
idates the use of this tool for the design of materials based on this genetic technique 12061 . Although 
the traditional approach is very useful in order to validate the results, genetic algorithms offer a way of 
designing crystal structures with desired properties, determining in this way the desired structure that 
underhes those properties. In an example where two types of colloids are grafted with complementary 
ssDNA sequences, ordered coUoidal crystals have been obtained. This methodology is generalizable, fast, 
and selective, able at the same time to reproduce parameters relevant to four currently realized crystals 
and unobserved structures 1206) . 


3.5. Atomistic simulations 

Atomistic simulations of 2 nm Au truncated octahedral NPs functionahzed with four single-stranded 
DNA chains have been also realized 1^12111 . However, such an approach is limited to a small number 
of short DNA chains similarly to the model of Starr and Sciortino 11421 . which deviates considerably 
from an experimentally realized nano-structure. However, in this study the authors have explored the 
possibility of hydrogen bonds between the adsorbed ssDNAs, finding that there is not a significant amount 
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of hydrogen bonds. This conclusion is consistent with the observed increase of the melting temperature 

(52l[m). 

Despite the extensive use of coarse-grained models, a million-atom MD simulation has been used to 
study the formation of bcc and fee super-crystals of DNA-functionalized Au NPs 12121 . The bcc crystal 
contained 2.77 million atoms, whereas the fee contained 5.05, with the diameter of the NP being 3 nm. 
After 10 ns DNA structures vary smoothly, while the ion distributions around DNAs can also be estimated. 
Analysis of these results shows that ions bind stronger in bcc structures than in fee structures. This effect 
explains the higher melting temperature of bcc structures and underlines the significance of entropic 
effects in fee structures. In this work, the Young’s and bulk moduli of the obtained super-crystals have 
been also obtained. Note that the Poisson ratios for both super-crystals were close to the ideal value, that 
is 1/3 (212). 

3.6. Simulations of anisotropic nanoparticles 

Anisotropic NPs grafted with DNA have recently appeared in the literature for the case of nano-cubes 
12131 . where the phase diagram similarly to the previous work with coated spherical NPs has been studied 
IM1I138I . These nano-cubes are capable of self-assembling into crystalline, liquid crystalline, rotator, or 
non-crystalline phases with both long-range positional and orientational order 12131 . 


4. Combined simulation and experimental studies 

The growth dynamics for DNA-mediated NP crystallization has been one of the fundamental ques¬ 
tions, because during the formation of crystalline-structure NPs become kinetically trapped. Moreover, 
it has been found that the growth of isolated crystallites is described by the power law ~ 12141 . The 

coalescence of these crystallites significantly slows down the growth process due to the orientational 
mismatch of the crystallites. MD simulations show that the misorientation angle decreases continuously 
during the coalescence, which is a signature of a rotation induced coalescence mechanism. The coarsen¬ 
ing of the conformation at the boundary between coalescing crystallites takes place after the orientation 
of the crystals. When nucleotide chains are larger, the dynamics become faster due to the enhanced sur¬ 
face diffusion, which more effectively reduces the curvature at the boundary of two super-lattices. These 
results have provided fundamental understanding on the relationship between NP surface chemistry and 
its crystalhte growth and coalescence 12141 . 

A quantitative prediction of the associated thermodynamics and kinetics of DNA-coated particles con¬ 
sidering different functionahzation schemes has been discussed by Leunissen et al. 12151 . Preparation, 
structural, and optical features of 2D cross-linked DNA/Au NP conjugates have been also considered in the 
literature I216II217I . as well as the tuning of ratios, densities and supra-molecular spacing in bi-functional 
DNA-modified Au NPs l218l . Experimental results (2151 show that the suspension behavior is very sen¬ 
sitive to the grafting details, such as the length and flexibility of the DNA strands, as well as the surface 
coverage. Therefore, in order to control the assembly process, emphasis should be put on the structural 
and dynamical features of the DNA coatings. By varying the temperature cycle, beads concentration, and 
surface coverage, the association-dissociation dynamics and kinetics can be studied systematically. More¬ 
over, the tethered backbone segments onto the NPs’ surface has been found to play an important role 
in the association-dissociation properties, since they affect the total number of inter-particle bonds and, 
as a result, the configurational entropy cost associated with these bonds. Independently of the tether 
backbone, self-complementary palindromic sticky ends can easily form intra-particle hairpins and loops 
considerably altering the association behavior. Hairpin DNA-functionalized Au colloids have found also 
application in the imaging of mRNA in living cells 12191 . Such structures are more important during faster 
temperature quenches, lower concentration and lower surface coverage. Other applications include the 
design of “self-protected” DNA-mediated interactions, which enable more versatile self-assembly schemes 

(ns). 

Self-replicating materials of DNA-functionalized colloids are also interesting 12201 . Such systems ex¬ 
ploit the strong specificity and thermal reversibility for the interactions between DNA-functionahzed NPs. 
The fundamental behavior of the self-replicating process related to the equilibrium and kinetic aspects 
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of aggregation-dissociation behavior of particles has been determined via experimental methods. It has 
also been found that the dissociation temperature is very sharp, occurring unexpectedly at low tempera¬ 
tures. When the fraction of the sticky ends becomes smaller, the dissociation temperature shifts to even 
lower scales. The sharpness of the transition and the fraction of sticky ends are important parameters 
in this self-replication scheme. Moreover, a double-stranded backbone prevents unwanted hybridization 
events, such as secondary structure formation, which leads to peculiar aggregation kinetics. This is due 
to the competition between inter- and intra-hybridization events. Additionally, by functionalizing a single 
particle species with two distinct DNA sequences, permanent bonds can form l220l . 

The bending dynamics of DNA-linked colloidal particle chains have been also studied 12211 by using 
a worm-like chain (WLC) model to represent the grafted DNA strands onto the NPs surface, analogously 
to the case of semi-flexible linear polymers. The persistence length of these grafted DNA chains can be 
determined by monitoring the thermal fluctuations of the chains. In this way, these lengths can be tuned 
within a large range from 1 to SO nm, which corresponds to a linking length of DNA of 7S to IS bases. 
Interestingly, the bending relaxation dynamics of these chains matches the theoretical predictions, sug¬ 
gesting the use of such a WLC DNA model approach for both equilibrium and dynamical studies 12211 . 

Park et al. 12221 investigated the structure of NPs aggregates at room temperature by combining theo¬ 
retical and experimental methods. Experimental methods involve extinction spectroscopy measurements 
and dynamic light scattering, where Au NPs have sizes between 60 and 80 nm and DNA strands of about 
30 base-pairs. Theoretical studies use calculated spectra for models of the aggregate structures to deter¬ 
mine the structure that matches the observations. Since particles in this study are larger than particles 
used previously, there is a higher sensitivity of spectra to aggregate structures. Combining theoretical 
and experimental results in the range where the highest agreement was found, it was argued that DNA 
hybridization takes place under irreversible conditions at room temperature. The model also takes into 
account other effects, such as the poly-dispersity in the size of NPs and lattice defects. The room temper¬ 
ature aggregate is always a fractal, and a morphological change from fractal towards compact structures 
occurs below the melting temperature 12221 . 

The kinetics of DNA-coated sticky particles through the exploration of a number of properties have 
been studied by Wu et al. (2231 . The aggregation rate, which depends on the rate of NPs encounter and 
the probability that such an encounter results in the sticking of the particle, has been studied. Contrary 
to the previous toy models (^, this one takes into account a possibility of hybridization. By combining 
experimental and theoretical techniques, the aggregation times of micron-sized particles as a function of 
DNA coverage and salt concentration have been obtained. A simple model using reaction-limited kinetics 
and experimental oligomer hybridization rates is in very good agreement with the experimental data. The 
Coulomb barrier at the nanometer scale retarding DNA hybridization is a significant controlling factor 
of the association process. This model also allows an easy measurement of microscopic hybridization 
rates from macroscopic aggregation, thus enabling the design of complex self-assembly schemes with a 
controlled kinetics 12231 . 


5. Synopsis 

To summarize, there has been a great advancement in understanding the self-assembly of DNA- 
functionalized NPs. A large number of crystalline structures can be fabricated and the effects of various 
parameters influencing the self-assembly are now well-understood. During the last years, a wealth of dif¬ 
ferent uses of DNA-functionalized NPs has been realized on the basis of different applications. However, 
the control of self-assembly and crystallization kinetics remains still an unresolved issue. This problem 
leads to kinetically trapped structures, one of the main reasons behind the difficulty in fabricating long- 
range stable crystalline structures and the efficient growth of crystalline structures on substrates (TtI . 
Problems in synthesis of these NPs in a controllable and predictable way is also a concern (iTl. Addition¬ 
ally, transitioning DNA-engineered NP super-lattices from solution to the solid state is challenging 12241 . 
Steps toward creating dynamic self-assembly materials, that is materials which are capable of chang¬ 
ing their crystalline structure in a controlled way due to specific stimuli, have been also taken 1T71I225I . 
Recently, dynamically interchangeable NP super-lattices have been produced via the use of nucleic acid 
allosteric effectors 12261 . This approach has been suggested for fee, bcc, CsCl and AIB 2 crystal symmetries 
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12261 . Another open question is the design of DNA-functionalized NPs that could lead to the formation of 
single-component non-compact diamond structures. Although some theoretical suggestions exist along 
these lines, this structure has not been realized experimentally, or by any simulation work. Moreover, 
given the large flexibility of DNA, the formation of single-component diamond structures with small NPs 
may not be possible 11561 . An opportunity window could be opened via the design with NPs of non- 
spherical shape which has been the topic of recent experimental I109II110I and simulation work l213l . 
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CaMOCKynMyeaHHJi DNA-(|)yHKi 4 i 0 Ha/ii 30 BaHi/ix ko/ioYaib 

n.E. TeoflopaKicP^, HJ. OiiacP, f. KayibP® H. /],e;i;iarcl® 

^ Biflfli/i xiMiHHOi iHxeHepiT, Ewnipia/i Koyieflx /lOHflOH, SW7 2AZ /lOHflOH, Be/iMKo6pMTaHm 
^ psHTp ripuK/iaflHMX MaTewaTi/iHHi/ix flOC/iiflxeHb, VHiBepci/iTeT KoBeHTpi, CV1 5FB, KoBeHipi, Be;ii/iKo6pi/iTaHifl 
^ iHCTI/ITyT Te0peTI/l4H0l c|)i3MKl/l, BifleHCbKHpi TeXHiHHMM yHiBepCl/ITei, A-1040 BifleHb, ABCTpia 
^ peHTp o6Hnc;iK)Ba/ibHoro Marepia/ioaHaBCTBae (CMS), A-1040 BifleHb, ABCrpia 
^ tpaxy/ibieT c|)i3MKi/i, YHiBepci/iTeT BiflHa, A-1090 BifleHb, ABCrpia 

Ko;ioi'AHi nacTMHKi/i 3 npi/iKpin/ieHi/iMi/i oahohutkobi/imh DNA(ssDNA)/iaHpK)>KKaMM, 3flaTHi caMOCKynnyBaTi/iCB y 
Hi/i3Ky pi3Hnx Kpi/icra/iiHHi/ix CTpyKTyp, fle ri6pnfli/i3aL4ifl /laHpioxKiB ssDNAcTBopioe 3B'fl3KM Mix KO/ioi'flaMi/i, CTa- 
6i;ii3yK)Hi/i Tx crpyKiypy. 3a;iexHO Bifl reoMerpiT i po3Mipy nacTi/iHOK, rycTHHi/i npi/ii±(en/ieHMX /lanpKDXKiB ssDNA, 
a Taxox rycTMHi/i ra BnOopy nocyiiflOBHOCTeii DNA, Moxna BuroTOByia™ Hi/i3Ky pi3Hnx Kpucra/iinHux crpyKiyp. 
OflHax 3a/iMLuaeTbra HeapoayMmi/iM i norpeOye iHieHCHBHi/ix flocyiiflxeHb ni/iTaHH^, akmm hmhom pi ((laKTopM 
poS/iaTb ci/iHepreTMHHuti BHecoxy npopec caMOCKynneHHfl DNA-cfjyHKpioHa/iisoBaHi/ix nacTHHOK HaHO- a6o Mi- 
Kpopo3MipiB. Bi/ibLiie Toro, flOflaiKOBi TpyflHOipi fl/iB po3B'B3aHHfl flaHoT npo6;ieMM no/iflraiOTb y BuroTOB/ieHHi 
fla/ieKOCflXHMX crpyKiyp 3aBflflKi/i KiHernHHUM kpmtmhhmm e/ieMeHiaM npw caMOCKynHeHHi. y crarri oOroBO- 
pKDKDTbCB HawHOBiiui flocnrHeHHfl TeopiT TB BKcnepuMeHTy 3 ocoOni/iBi/iM Haro/iocoM Ha ocraHHi floc/iiflxeHHfl 
Meroflawi/i KOMn'roiepHoro MOfleuKDBaHHB. 

KnioHOBi cnoBa: DNA-(j>yHKu,i0Hani30BaHi HaHO>^acTi/iHKn, caMOCKyn^enHa, eKcnept/iMenT, reopifi, 
KOMn'miepne MOflejitOBaHHa 
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